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ABSTRACT 

We examine the status of the threat posed to stellar-mass black hole candidates by the 
possible existence of Q-stars (compact objects with an exotic equation of state which 
might have masses well above the normally-accepted maximum for standard neutron 
stars). We point out that Q-stars could be extremely compact (with radii less than 1.5 
times the corresponding Schwarzschild radius) making it quite difficult to determine 
observationally that a given object is a black hole rather than a Q-star, unless there 
is direct evidence for the absence of a solid surface. On the other hand, in order for 
a Q-star to have a mass as high as that inferred for the widely-favoured black hole 
candidate V404 Cygni, it would be necessary for the Q-matter equation of state to 
apply already at densities an order of magnitude below that of nuclear matter and 
this might well be considered implausible on physical grounds. We also describe how 
rotation affects the situation and discuss the prospects for determining observationally 
that black hole candidates are not Q-stars. 
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1 INTRODUCTION 



The search for an unequivocal demonstration of the exis- 
' tence of stellar-mass black holes has focussed on X-ray emit- 
, ting binary systems consisting of an ordinary star together 
■ with a compact object. If the mass of the compact object (as 
I determined from kinematical measurements) appears to be 
greater than the maximum possible for a neutron star Mmax, 
then the object has been acknowledged as a black hole candi- 
date. The value of Mmax is still not reliably known, because 
of uncertainties in the equation of state of matter at high 
densities, but it has been widely believed that the limit of 
3.2 Mq derived by Rhoades & RufRni (1974) (together with 
a possible 25% upward correction for rotation) gives a secure 
upper bound. 

The best current stellar-mass black hole candidates are 
in soft X-ray transients (SXTs), a sub-class of the low mass 
X-ray binaries (see Charles 1996). During quiescence, the 
accretion disc in these systems becomes extremely faint and 
it is then possible to carry out detailed photometry and 
spectroscopy of the optical companion, which allows the 
mass of the compact object to be directly determined (van 
Paradijs & McClintock 1995). Fig. 1 shows the presently- 
known masses of neutron stars and black-hole candidates. 
The "neutron star" masses all lie within a small range of 
1.4 Mq, whereas the black-hole candidates seem to form 
a distinct grouping around ~ 10 Mq. The most convinc- 



ing of the black hole candidates are V404 Cyg (Shahbaz 
et al. 1994b) and Nova Sco (Orosz & Bailyn 1996) with 
masses of 12 Mq and 7 Mq respectively, well above the 
Rhoades/Ruffini limit. 

Rhoades and RufBni derived their result in response to 
the problem caused by different high-density equations of 
state leading to widely different values for Mmax, the idea 
being to derive a firm upper limit for non-rotating mod- 
els on the basis only of knowledge which could be consid- 
ered as completely secure. However, some of the assump- 
tions made are, in fact, distinctly questionable (see Har- 
tle 1978, Friedman & Ipser 1987). In particular: (i) it was 
assumed that the equation of state can be taken as accu- 
rately known for densities up to a fiducial value po which 
they took as 4.6 x 10^^ gcm~^; (ii) they imposed a causal- 
ity condition which would only be appropriate for a non- 
dispersive medium. If a more conservative value is taken for 
po (10^* g cm^'^), the causality condition is dropped and al- 
lowance is made for rotation, the upper bound for Mmax 
obtained in this way goes up to 14.3 Mq (Friedman & Ipser 
1987) which is no longer useful in considering black hole 
candidates such as V404 Cyg. However, standard realistic 
neutron star equations of state do, in practice, give masses 
satisfying the original condition A4 < 3.2 Mq even with ro- 
tation. (It should be noted that all of this discussion of the 
mass limit is being made within the context of general rel- 
ativity which is taken to provide a correct description of 
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gravity within this strong-field regime. We work within this 
context throughout the present paper. However, while gen- 
eral relativity is well-accepted as our theory of gravity, we 
recall that it has not yet received convincing experimental 
verification away from the weak-field limit and so one should 
not overlook the possibility that this might not be the cor- 
rect description.) 

How much can the standard neutron star equations of 
state be trusted for matter at high densities? Is it correct 
that neutron star matter consists of a mixture of neutrons, 
protons, electrons, mesons, hyperons, etc., and is held to- 
gether just by its own self-gravity? Quantum chromodynam- 
ics (QCD) contains the idea of confinement by the strong 
force which is normally thought of in terms of quarks be- 
ing confined within nucleons. It was in connection with this 
that the strange star model was introduced (Witten 1984; 
Haensel, Zdunik & Schaeffer 1986; Alcock, Farhi & Olinto 
1986), resting on the hypothesis that strange quark mat- 
ter might bo the absolute ground state of baryonic matter 
oven at zero pressure. Strange stars would be essentially sin- 
gle giant nucleons with baryon number A ~ lO^'^ and with 
the confined quarks being free to move within the false vac- 
uum which extends throughout the interior. They could have 
masses and radii similar to those of standard neutron stars 
and have been advocated as a viable alternative model for 
pulsars although there is a problem over explaining glitches. 

Strange stars are "safe" as far as the mass limit is con- 
cerned. Although their equation of state is very different 
from that for standard neutron star matter, the maximum 
mass is within the standard range (it is ~ 2.0 M© for non- 
rotating models with some variation depending on uncer- 
tain parameter values). However, some effective field theo- 
ries of the strong force allow for it not only to confine quarks 
in the normal way but also to confine nucleons (neutrons 
and protons) at densities well below that of nuclear matter 
[pnm ~ 2.7 X 10"'^'* gcm~^), giving an equation of state dif- 
ferent from the standard one at densities below the values 
normally taken for po. Models based on this idea were in- 
troduced by Bahcall, Lynn & Selipsky (1989a,b, 1990), who 
named them Q-stars (although note that the "Q" here does 
not stand for "quark" but for a conserved particle number) . 
These are not safe for the mass limit and might, in principle, 
have very high masses up to more than 100 M©. 

Even if one rules out Q-stars as models for pulsars (they 
have similar difficulties in this respect as for strange stars) 
there still remains the possibility that they could be an al- 
ternative model for the more massive objects in black-hole- 
candidate systems. In connection with this, there are a num- 
ber of questions which immediately present themselves. For 
a given mass, how small could a Q-star be? How physically 
reasonable are the versions of the Q-star equation of state 
which would allow masses as high as that of the compact 
object in V404 Cyg? What happens when rotation of the 
object is considered? These issues are considered in the next 
Section. 

2 PROPERTIES OF Q-STAR MODELS 

For our calculations, we have used the simplest form of the 
Q-star equation of state: 



where p is the density, p is the pressure, f/o is the energy 
density of the confining scalar field and Uv measures the 
strength of the repulsive interaction between nucleons. This 
represents chiral Q-matter (for which the particles have zero 
mass within the false vacuum) but the results obtained are 
only marginally different for non-chiral Q-matter. It is con- 
venient to introduce the parameter (= avU^^'^ it / since, 
when this is fixed, all results scale with Uq. 

Fig. 2 shows representative mass/radius curves for Q- 
stars and for a standard neutron star equation of state 
(model C from the collection of Arnett & Bowers 1977) 
and it is clear that the two curves are quite different. The 
lowest-mass Q-star models, for any given C,, have an almost 
constant value of the density, giving M oc _R'\ but with 
increasing mass, the density profile becomes progressively 
more peaJsed towards the centre. As the central density is 
increased further, a maximum mass is eventually reached 
and this marks the end of models which are stable under ra- 
dial perturbations. The maximum-mass model is the most 
compact stable one (i.e. R/M is a minimum) and therefore 
the closest in size to a black hole of the same mass. In Fig. 3, 
we show {R/M)min plotted as a function of C for Q-stars to- 
gether with corresponding values for a black hole and for the 
representative neutron star equation of state. (We are us- 
ing the standard geometrical units of general relativity with 
c = G = 1.) It can be seen that Q-stars can be very com- 
pact with R/M < 3 (although note that very large values 
of are not physically reasonable.) A stable, non-rotating 
12 Mq Q-star might have a radius as small as ~ 52 km as 
compared with ~ 36 km for a non-rotating black hole of 
the same mass. Having R/M < 3 means that the surface 
lies inside the circular photon orbit, the place where cen- 
trifugal force becomes attractive rather than repulsive (see 
Abramowicz & Prasanna 1990), but the models do not seem 
to be quite compact enough to allow for some more exotic 
relativistic effects such as resonance of axial gravitational 
wave modes (Chandrasckhar & Ferrari 1991) or production 
of an internal ergoregion when the object is set into rapid 
rotation (see Butterworth & Ipser 1976), which would lead 
to an instability. 

The maximum mass which a Q-star could have is related 
to the threshold minimum density, Pmin, above which the 
Q-matter equation of state is taken to apply. (This is the 
density at which the pressure goes to zero.) The relation is 
shown in Fig. 4 for the limiting cases C = ^-^d oo 
and it can be seen that varying C, over the entire range in 
between makes only a small difference. In order to have a 
non-rotating Q-star with as high a mass as that inferred for 
the compact object in V404 Cyg (marked with the horizontal 
dashed line) , it would be necessary for the threshold density 
to be below p„rn (marked with the vertical dashed line) by 
about a factor of 10 or more. Rotation makes little difference 
to this. Even if one admits the idea of Q-matter in principle 
for some range of densities around pnm, such a very low 
value of Pmin might seem implausible. 

Compact objects accreting matter in binary systems 
will be spun up by the angular momentum of the accreted 
material and so it is important to ask what effect this would 
have for the comparison between the black hole and Q-star 
pictures. Fig. 5 shows the locations of the equatorial radii 
of a black hole, an extreme Q-star (^ — » oo, maximum com- 
pactness) and a representative neutron star, as functions of 
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a/M whore a is the angular momentum per unit mass and all 
quantities are again in geometrical units. All of the radii are 
calculated in the same coordinate-independent way, by tak- 
ing the proper distance round the circumference of the circle 
concerned and dividing by 2n. Doing this, the equator of the 
black hole event horizon is always at r = 2M irrespective 
of the value of a/M. (We emphasize that this is in contrast 
with the situation for the commonly-used measure in Boyer- 
Lindquist coordinates for which the equator is at r = 2M 
when a = but r M as a ^ M. Other frequently-quoted 
quantities such as observed rotational frequencies and bind- 
ing energies are usually calculated in an invariant way.) The 
formulae used for the Q-star and neutron star models arc 
calculated within the slow-rotation approximation, correct 
to second order in the rotational velocity Q and comparing 
rotating and non-rotating models with equal central den- 
sity. Also shown in Fig. 5 is the location of the marginally 
stable orbit (corresponding to the inner edge of a Keplerian 
accretion disc) which, for a given a/M, is the same for all of 
the objects within the slow rotation approximation. While 
bearing in mind that results obtained with the slow-rotation 
approximation should only be regarded as indicative for the 
higher values of a/M, several interesting conclusions can be 
drawn from this figure. As a/M is increased, the r/M of the 
marginally stable orbit decreases and eventually it reaches 
the equator of the neutron star which expands out to meet 
it. For the extreme Q-star, r/M of the equator actually de- 
creases slightly and the marginally stable orbit would reach 
it only for larger a/M. 



3 DISCUSSION AND CONCLUSIONS 

How can one be sure that a high-mass compact object such 
as that in V404 Cygni is a black hole and not a Q-star? Any 
evidence for the existence of a solid surface {e.g. by observ- 
ing a Type I X-ray burst) would, of course, immediately rule 
out the possibility of the object being a black hole but in 
the absence of such evidence, what can be done? Narayan, 
McClintock & Yi (1996) have proposed a model for some 
SXTs in quiescence which is in agreement with available 
observational data and depends on the compact object in 
question being a black hole. However, this is still an indirect 
argument. What one would really like is to have direct obser- 
vational evidence of accreting material at radii smaller than 
would be possible with a Q-star. With satellite experiments 
such as the Rossi X-ray Timing Experiment (RXTE) and 
the Unconventional Stellar Aspect experiment (USA) being 
capable of time resolution down to the order of 1 /isec, there 
is a possibility of detecting the location of the inner edge 
of the accretion disc (if, indeed, there is a well-defined in- 
ner edge) by seeing a corresponding frequency cut-off in the 
power spectrum, by seeing dying pulse trains from bright 
features as they reach the inner edge and fall in (Stoeger 
1980) or from analysis of quasi-periodic oscillations (Miller, 
Lamb & Psaltis 1997; see Zhang et al. 1997 for discussion 
of a possible detection of the inner edge of the disc by this 
means for binaries with compact components in the neu- 
tron star mass range). For non- rotating compact objects, 
the marginally stable orbit is at r = 6M, much larger than 
the radius of the most interesting Q-star models, but with 
increasing rotation of the compact object, it moves inwards 



and if the compact object were a Q-star, it would even- 
tually meet the surface. If evidence were to be found for 
an accretion flow extending further inwards than would be 
possible with a Q-star, this would then point very strongly 
towards the compact object being a black hole, since the 
Q-star model probably presents the last realistic possibility 
for avoiding that conclusion within the context of general 
relativity. 
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FIGURE CAPTIONS 

Figure 1: Mass distribution of neutron stars and black 
holes for which masses have been directly measured (Gies 
& Bolton 1986; Nagase 1989; Shahbaz, Naylor & Charles 
1993; Thorsett et al. 1993; Shahbaz, Naylor & Charles 1994a; 
Shahbaz ct al. 1994b; van Paradijs & McChntock 1995; 
Beekman et al. 1996; Orosz & Bailyn 1996; Beekman et al. 
1997; Remillard et al. 1996; Reynolds et al. 1997; Shahbaz, 
Naylor & Charles 1997). The components of the binary pul- 
sars are marked with with P (pulsar) and C (companion). 
Also shown are the X-ray pulsars. The neutron star systems 
lie at ~ 1.4 Mq, whereas the black- hole candidates seem to 
cluster around ~ 10 Mq. The strongest black hole candi- 
dates are V404 Cyg and J1655-40. The vertical hashed line 
represents the Rhoades/Rufflni limit of 3.2 Mq. 

Figure 2: Representative mass/radius relations for (a) Q- 
stars and (b) standard neutron stars. Note that the scaling 
and precise form of the Q-star curve depend on parameter 
values which uncertain. 

Figure 3: {R/M)min, the minimum value of (R/M) for sta- 
ble non-rotating Q-star models, is plotted as a function of C- 
Also shown, for comparison, are values for a representative 
neutron star equation of state and for a black hole. For the 
Q-star, {R/M)min tends towards a constant value of 2.8 as 
C — » oo: Q-stars can be very compact. 

Figure 4: The maximum mass of non-rotating Q-star mod- 
els is plotted as a function of the threshold density for Q- 
mattcr pmin- Curves are drawn for the limiting cases C = 
and ^ —^ oo. The horizontal dashed line corresponds to the 
mass of the compact object in V404 Cyg; the vertical dashed 
line marks the rmclcar matter density pnm- To allow for a 
Q-star with a mass as high as 12 Mq, the threshold density 
would have to be about an order of magnitude below nuclear 
matter density, or less. 

Figure 5: Equatorial radii of a black hole, an extreme Q- 
star, a typical compact neutron star and the location of 
the marginally stable orbit are plotted as functions of a/M. 
These are (coordinate independent) circumferential proper- 
distance radii calculated consistently within the slow rota- 
tion approximation. The equatorial radius Req of the black 
hole measured in this way is constant, independent of a/M. 
Req/M for the most extreme Q-stars actually decreases very 
slightly as a/M is increased within the slow rotation regime, 
which is in contrast with the behaviour of the neutron star. 
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